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Q\ Abstract 



A perturbative technique, the low-temperature expansion, is de- 
veloped for matrix models of random surfaces. It can be applied to 
models with arbitrary target spaces, including ones with c > 1. As a 
^5 simple illustration, the series is worked out to 10th order for the sur- 

lf} face coupled to a g-state Potts model. Accurate estimates for, e.g., j s t r 

^ — I are obtained both in the low q (c < 1) and high q (branched polymer) 

regimes, including the logarithmic corrections to scaling. 
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^ 1 Introduction 

in 

The trouble with matrix models of 2D random surfaces and quantum gravity 
[1] is that "we are in the unfortunate situation of either solving exactly (for 
+ 7 > c < 1) or not understanding at all what is happening" [2]. If the target 

Qh space has loops, we do not even know where to begin. 

In this paper I develop a perturbative technique, the low temperature 
expansion (lte), which can handle any multi-matrix model, and therefore 
any kind of target space or matter model on the surface [3, 4]. In effect, 
the LTE expresses multi-matrix model correlation functions in terms of one- 
matrix model correlation functions. The result is a series in powers of e _/3 
((3 is the inverse matter temperature), whose coefficients depend on the 
cosmological constant, g, and on the geometry of the target space (which I 
will always take to be a graph on a discrete set of points). In fact, most of 
the calculations required to obtain these coefficients are universal — they are 
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identical for any target space. The target space enters at only one step, in 
the calculation of so-called "homomorphism coloring factors." 

There exists another series technique for matrix models, the strong- 
coupling expansion, where the expansion parameter is the cosmological con- 
stant [5, 6, 7]. There, the n-th term includes all surfaces made up of n 
simplices. The difficulty of the strong-coupling approach — which has been 
pushed to about 10th or 20th order, depending on the target space — is that 
one cannot always trust extrapolations from such small surfaces, especially 
above c = 1, where the geometry becomes quite singular. 

The LTE does not suffer from this difficulty. In every order it includes 
surfaces of all sizes. The lte is a spin cluster expansion: what is constrained 
is not the total area of the surface, but the number of spin clusters. 1 With 
very moderate computational effort it gives accurate results for the critical 
points and exponents (including the logarithmic corrections) for c < 1, as 
well as for the (even more) branched large-c surfaces. 

In the present article, following a detailed discussion of how to calculate 
the LTE 2 , I carry the series out to 10th order for the simplest multi- matrix 
model, the g-state Potts model. 3 The n-th term of the lte, the coefficient 
of e~ n/3 , is now an explicit function of g and q. In other words, one does not 
have to do a new calculation for each new value of q — once the calculation of 
the LTE is finished, one simply plugs in values of q into the series coefficients. 
The Potts example is meant foremost as an illustration of the effectiveness 
of the LTE; but also as a study of a model interesting in its own right. What 
follows is a brief discussion of the Potts model on a random surface: first of 
all, why is it interesting? 

A standard argument for why it's not runs as follows: to get a non-trivial 
modification of pure gravity, one must take a model which — on a fixed 2d 
lattice — has a continuous transition, and couple it to the surface. But for 
q > 4, the Potts model has only discontinuous transitions on a fixed lattice. 
Doesn't this mean that putting a q > 4 model on a random surface should 
yield only pure gravity? No, because coupling a spin model to a random 
surface (usually) softens the spin model's phase transition. For instance, 

x The number of links between spin clusters, to be exact; but the number of clusters is 
never greater than one plus the number of links between them. 

2 The lte was introduced in [3], but as a qualitative rather than quantitative method. 
In [4, 8] it was resummed to all orders in the q, c, d — > oo limit. 

3 "Simplest" because its target space, K q , has no structure whatsoever: every point (or 
matter state) is connected to every other point, with the same weight. Consequently, the 
homomorphism coloring factors reduce to ordinary chromatic polynomials — see section 2.5. 
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the transition of the q = 2 Potts model (the Ising model) is second order 
(a = 0) on a fixed 2d lattice, and third order (a = — 1) on a random surface. 
The effect is even more dramatic at q = oo: on a fixed lattice the model is 
very first-order, while on a random surface it is third-oidei. At q = oo the 
model is in fact exactly and explicitly solvable [8]; surprisingly, it is identical 
to c = oo models (such as the d = oo lattice, or infinitely many Ising spins) 
[4]. This, then, is the interest of the Potts model on a random surface: for 
q < 4 it corresponds to c < 1, while for q = oo it is identical to c = oo. 
Can the 4 < q < oo region teach us anything about c > 1? And if there are 
continuous phase transitions for q > 4 that modify pure gravity, do these 
correspond to some new flat models coupled to gravity, or are these models 
unique to random surfaces? 

The model is defined as follows. There are two kinds of degrees of free- 
dom: the lattice and the matter. The lattices are randomly triangulated 2D 
spherical surfaces, which are dual to spherical cubic fat graphs. The matter 
consists of a spin a a 6 {1, . . . , q} on each triangle a (or on each vertex of 
the fatgraph). For a given geometry, the matter Hamiltonian is 



where the sum runs over all edges of the fatgraph (a triangle may be con- 
nected to itself, and two triangles may be connected by more than one edge). 
The partition function of the model is realized by the following matrix inte- 
gral (putting e~ 2/3 = z): 



where the fa are N x N hermitian matrices and the matrix integral is to 
be understood as being divided by N 2 and taken in the limit N — > oo. 
The matrix T defines the matter that is to be coupled to the fluctuating 
surface — it is the connection matrix of the target space graph; in the case 
of the Potts model this is 



It will be convenient to use a rescaled coupling constant a = Az, A = q — 1. 
Finally, the normalization f is defined by requiring F q (0,z) = 1. 



(i) 





(3) 
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What do we expect to find? First, for every value of q there should be a 
curve of critical points in the (g, a) plane. For q < 4 (c < 1) we know that 
this curve has two phases: the low-temperature magnetized phase and the 
high-temperature disordered phase, both of which are in the same univer- 
sality class as pure gravity (i.e., they both have ^ str = — |). The boundary 
between these phases is the critical point of the spins, where their critical 
fluctuations modify the global geometry, and at which 7 sir / — |, in gen- 
eral. For very large q the picture is different: there is still a low-temperature 
phase (where ^y str = — |), but this now goes into a high-temperature phase of 
branched polymers (where j str = +|); at the multicritical point separating 
the two phase ^y str = +|. As q is decreased, the old disordered phase makes 
its re-appearance 4 ; the order of the phases is: magnetized at low temper- 
ature; branched polymers at intermediate temperature; and disordered at 
high temperature. As q decreases the magnetized and the disordered phase 
move closer and closer to each other until, at some q = q c , they touch: 
the branched polymer phase is eliminated, and the low-g phase diagram is 
restored [9]. 

Here is a summary of the results that I get using various series extrap- 
olation techniques on the 10th order Potts model series. First, I calcu- 
late, for different values of q, 7 sir on the critical (</, a) curve using Pade 
approximants — this is plotted as a function of g in Fig. 4. For small q there 
are two ^y str = — | regions, the magnetized and the disordered, separated by 
a sharp peak at the phase transition (the height of the peak agrees quite 
well with KPZ for q = 2 and 3). As q gets larger, the peak gets broader (an 
effect also observed in simulations). For high enough q, moreover, a plateau 
seems to appear between the two regions, its level approaching j str = +| in 
the limit q — > oo. To accurately locate the phase transition for q > 4 (where 
it may no longer correspond to the j str peak), and to include logarithmic 
corrections, I have used a variation on the ratio method. The results are 
summarized in Fig. 5 and in Table 3. The logarithmic corrections give a 
spectacular improvement in 7 sir at and just above q = 4. For example, for 
q = 4 (where KPZ gives 7 str = 0), I get 7 sir pz —0.2 without logarithmic cor- 
rections, but j str « 0.003 with logarithms. There is a region, approximately 
4 < q < 15, where j str increases gradually, while the logarithmic exponent 
remains roughly constant. Then, for higher q, the logarithmic exponent 
and the corresponding 7 str decrease, probably indicating that logarithmic 
extrapolations can no longer be trusted in this regime (where the logarith- 

4 In the language of [9] this is called the phase of "large galaxies." 
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mic corrections are small or zero). Above this uncertain region (roughly 
15 < q < 50), it seems clear that one should use non- logarithmic fits. The 
resulting j str increases slowly and smoothly with q, approaching 0.37 in the 
limit q — > oo (the exact answer is j str = +|). 

The rest of this paper is organized as follows. Section 2 is a careful 
discussion of the LTE (for arbitrary target spaces) and how to calculate 
the coefficients. Section 3 is a numerical analysis of the Potts example. 
Conclusions and prospects for future work are all to be found in Section 4. 

2 Generating the series 

This section is a didactic presentation of the mechanics of the lte. I have 
tried to give motivations and intuitive explanations of the different factors 
that enter. Table 1 summarizes the calculations though fourth order. There 
are two types of factors: ones that do not depend on the matter model that 
is coupled to the surface (sections 2.1-2.4; all but the right-hand column of 
Table 1); and one that does (section 2.5; the right-hand column of Table 1). 

2.1 Spin clusters 

The low-temperature expansion for discretized random surfaces is similar to 
the LTE for ordinary spin models. It is essentially a spin cluster expansion: in 
each order a n (recall that the matter coupling constant a = (q — l)e~ 2/3 — see 
eqs. (2) and (3)) one sums all configurations that have n links that connect 
unequal spins. It is convenient to represent these configurations by skeleton 
graphs (or "skeletons" for short). The vertices of a skeleton are the spin 
clusters, that is, connected sets of lattice sites that have a constant spin; 
while the edges of a skeleton are only those edges of the lattice that connect 
the spin clusters to each other. The spin clusters are surfaces of arbitrary size 
(whose fluctuations will be summed by appropriate vertex factors), but with 
restricted topology. (The reader should not confuse the discretized random 
surfaces — i.e., the fatgraphs of a matrix model — with the skeleton graphs 
of the lte. A skeleton represents a fatgraph with a particular configuration 
of spins; many such fatgraph/spin configurations can be represented by one 
skeleton.) By definition, a vertex of a skeleton cannot be connected to itself. 

For instance, the skeletons that enter into orders 0-4 of the LTE are given 
in the left column of Table 1. Weighted by different embedding factors, 
these graphs should occur in the LTE of any model — random surface or fixed 
lattice. For any model on a fixed, 2d square lattice however, only two of 
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Table 1: Graphs through order 4 of the LTE and the factors that they con- 
tribute 
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the graphs in Table 1 occur: • and ' >T (for example, • — • does not 
occur because it is impossible to divide the lattice into two clusters that 
are connected by just one edge). In fixed-lattice models, each skeleton typi- 
cally gives rise to three factors: symmetry, embedding, and "homomorphism 
coloring" (the number of ways to assign spins to the clusters). 

For a random surface model, almost all of the graphs in Table 1 will 
typically contribute. Each skeleton contributes two of the above factors: 
symmetry and coloring. Now, however, we have a random geometry model. 
Therefore the spin clusters — the vertices of the skeleton graph — are fluctuat- 
ing random surfaces, and the contribution from each skeleton must include 
sums over these fluctuations. These will appear as vertex factors, denoted 
as 7r's or p's in what follows. Each cluster will be a punctured surface, with 
as many punctures as the degree of the corresponding skeleton graph vertex. 
Since the spins on each cluster are frozen, the clusters' geometrical fluctua- 
tions will be governed by a pure gravity one-matrix model. Finally, there is 
no restriction on the area of the spin clusters: each order of the lte includes 
surfaces of all sizes. 

Notice that only the coloring factors for the skeletons depend on the 
details of the matter model that is coupled to the surface. I will return to 
the coloring factors in section 2.5, while the next sections will be devoted to 
the topology and geometry of the spin clusters. 

2.2 Twist and topology 

Exactly which surfaces must be summed over for each spin cluster? The 
answer depends, unfortunately, on the details of each skeleton graph. I will 
first introduce a classification for punctured discretized random surfaces, and 
then explain which of these contribute to any cluster in any given skeleton. 

Closed discretized random surfaces are classified by their genus, or the 
minimum number of handles. The additional complication will be that spin 
clusters, since they are (usually) connected to other spin clusters, are open, 
punctured, surfaces. In the context of the lte, a natural way to classify 
punctured surfaces will make use of two parameters: the genus and the 
twist. The genus is defined by putting endcaps on the punctures, so that 
none of the endcaps is connected to any other (see Fig. la). The genus of 
the punctured surface will be defined as the genus of the endcapped version. 
Our clusters will all be genus-0 surfaces according to this definition. 

To define the twist, put a dummy matrix A on each puncture. The 
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(a) Definition of genus 



(b) Definition of twist 



TrA 



TrA TrA 



(TrA) 



Figure 1: Classification of punctured surfaces 




(a) 



(b) 



Figure 2: An example of (a) an untwisted spin cluster (counted by 7r 2 ) and 
(b) a twisted spin cluster (counted by tt^i) with two punctures 
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Figure 3: A planar surface with a twisted cluster. The skeleton graph is 
shown on the left, an example of a possible corresponding surface (fatgraph) 
on the right. 

propagators will connect the A's in some way; the result of these contractions 

will be (see Fig. 1) Tr A kl Tr A* 2 • • •, so that k 1 + k 2 -\ = k, where k is the 

number of punctures; the set {ki, k 2 , . . .} is therefore an (integer) partition 
of k. 5 For a surface with k punctures, the quantity which I will call "twist" 
will be the corresponding partition p chosen from P k , the set of partitions of 
k; for example, P 3 = {{3}, {2, 1}, {1, 1, 1}} (see Fig. lb). Fig. 2 shows the 
simplest concrete example of twist, for k = 2. I will denote the corresponding- 
Green's functions by 7r p (g), so that the three possibilities for k = 3 in Fig. 
lb will be called (from top to bottom) 7r 3 (</), 7r 2i i(g), and 7r liljl (g). 

For a cluster with k punctures, call the surfaces counted by n k untwisted, 
and all the others (counted by Kk-i,i, ■ • •) twisted. The untwisted Green's 
functions n k are usually called just the "Green's functions"; the twisted 
surfaces, to my knowledge, have never been considered or calculated (see 
section 2.3 for how it's done). Your intuition tells you that twisted surfaces 
cannot contribute to the planar limit; this intuition is false. To properly 
represent multi-matrix models in terms of clusters, we must include some 
twisted clusters. The reason is that, in some cases, twists can be untwisted: 
the resulting surface is really planar, and therefore must be included. Fig. 
3 provides a simple illustration. 

Now we can return to question of which skeleton graphs contribute to 
the partition function. Since we want to have a planar surface in the end, 
the skeleton itself should be thought of as a fatgraph. Each vertex (spin 
cluster) with k legs contributes a factor of k\~ 1 n p (g), where p is a twist, 
i.e., a partition of k: this vertex factor is the contribution of the fluctuating 
geometry of the spin cluster. Thus, k-\eg vertices come in w k flavors, where 
w k is the number of integer partitions of k (for large k, w k grows as e^). In 
principle, each &-leg vertex contributes a vertex factor which is a sum over 
all the twists of k with equal weight, ^2 pePk 7r p . However, some combinations 
of twists will result in non-spherical surfaces; so that the net vertex factor 

6 An (integer) partition of k is an orderless set of positive integers whose sum is k. 
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is actually some linear combination of the twists, J2 P eP k c p n p — see Table 1, 
column "Spin cluster /topology" for examples. 

In calculating the linear combinations of twists that occur at each vertex 
of each skeleton, untwisted clusters n k may be faithfully represented as ordi- 
nary linear matrix model vertices (Tr A k ), while the twisted clusters n kltk2t ... 
are represented as nonlinear vertices (popularized by [10]) Tr A* 1 Tr A* 2 • • • 
which appear in models of polymerized surfaces. In other words, a twisted 
spin cluster is topologically equivalent to a surface touching. This is why the 
nonlinear vertices occur in the effective potentials of the one-matrix version 
of the Potts model [11], as well as in the effective potential of the reduced 
lattice model [12]. 

The preceding calculation may be simplified by noting that only certain 
linear combinations of the 7r's occur as vertex factors. To determine which 
linear combination to take for a given vertex of a given fatgraph, let us 
define the "connectivity partition" of a vertex. Consider a vertex with k 
legs in some skeleton. Any two of its legs are said to be connected if one 
can reach one of the legs starting from the other, without going through the 
vertex itself. The k legs will thus be divided into some number of clusters: 
every leg is connected to every leg in its cluster, and not connected to any 
other leg. If the number of legs in the clusters is qi, q 2 , . . ., then the q's are a 
partition of k, the connectivity partition. If the graph is tree, for example, for 
a given vertex no leg will be connected (in the present sense) to any other, 
and therefore the connectivity partition of every vertex will be {1, 1, . . .}. 
This connectivity partition precisely determines the linear combination of 
7r's which will be the required vertex factor. I will denote these vertex factors 
p qi ,q 2 ,...- Here are some examples: 
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, ^2,2 

P2,2 = 7I"4 H Y~ 

, ft 3,1 , n 2,2 , 7r 2,l,l 

P2,i,i = *, + — + — + _ 

Pi, 1,1,1 = ^4 + 7T3.1 + 7r 2 , 2 + 7I"2, 1,1 + ^1, 1,1,1 

The recipe for using the p vertex factors is as follows. Instead of summing 
over all ^-twists for each A;- vertex, take only the untwisted vertex, Trtp k . 
Having calculated the large- N factor using the untwisted vertices, determine 
the connectivity partition q of each vertex, and multiply by p q . Using the 
p rather than ir vertex factors makes the calculation more transparent and 
saves one the trouble of summing over all twists for each vertex. 



2.3 Calculating the clusters 

The spin cluster Green's functions n p are generated by an external matrix 
integral [11, 13] 

n(A) = - log J V<p exp -NHt [Kip + V(tp)] (5) 

where A is some fixed N x N hermitian matrix, and as usual matrix integrals 
are understood in the sense lim^r^oo N~ 2 (and tv X = TrX/ dimX). The 
K(p vertex inserts punctures, and the fact that A is a matrix keeps track of 
the twist. Expand the generating function 11(A) in powers of A, one finds 
precisely 

11(A) = 7T + 7Titr A + ^[7r 2 tr A 2 + 7r M (tr A) 2 ] 

+^[7r 3 tr A 3 + 7r 2 ,itr A 2 tr A + 7r 111 (tr A) 3 ] + • • • (6) 

For the cubic potential V(ip) = cp 2 /2 + gcp 3 the generating function 11(A) 
has been calculated exactly [11, 13]. One could expand this result in powers 
of A to obtain closed-form expressions for the -k p as functions of g. 

A more efficient method is to use the Schwinger-Dyson equation for 
integral (5), which reads 

A-Nd A U + 3g[d 2 A U + N 2 (d A U) 2 } = (7) 

Inserting the expansion (6) into this equation, one finds in n-th order a 
system of linear equations for the n-puncture functions n p (p 6 P n ), where 
the coefficients and inhomogeneous terms involve functions with fewer than 
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n punctures; this gives us, therefore, a recursive procedure for calculating 
the 7T P . The only complication is that the most twisted function in each 
order (7Ti ... i) does not appear. This is not a problem, though: one can 
easily calculate the sum over all twists in each order Tr n = ^2 pePn tt p by 
setting A = Al; eq. (6) now reads 

n(A) = 7r + 7r 1 A + ^A 2 + ^A 3 + --- (8) 

by simply shifting the matrix in the BIPZ one-matrix integral <p — >■ <p + xl 
by a multiple of unity; this gives us an extra linear equation for 7Ti ...i. 

The Schwinger-Dyson equations (7) are valid for all N, and therefore 
generate non-spherical corrections; for the sake of efficiency, it would be 
desirable to eliminate these corrections. This can be done by following a 
certain prescription for the matrix derivative operators d\ and d\. For a 
partition p = {p 1 , . . . define \p\ = k; p—pi as p with the l-ih element re- 
moved; and x P (A) = A Pl • • • A Pk . The prescription for the matrix derivatives 
is 

bl 

SaX p (A) = J^Pr^Xp-pM) (9) 

r=l 

bl /p-- 2 \ 
d 2 A x P (A) = J2Pr k p r- s - 2 Tr A s X P - Pr (^) + non-spherical (10) 

The first derivative (9) is exact; the non-spherical terms omitted in the sec- 
ond derivative (10) are guaranteed to contribute only subleading 0(N~ r ) 
terms to the Schwinger-Dyson equations; moreover, all the subleading terms 
in the Schwinger-Dyson equations are eliminated by simply omitting the 
"non-spherical" terms in eq. (10). Putting eqs. (7), (6), (9) and (10) to- 
gether, the first few equations of the Schwinger-Dyson hierarchy are 
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3<7(7T 2 + 7T?) - 7T! = (11) 

j + 27ri7r 2 ) - tt 2 + 1 = (12) 
(? + T +27ri7ri - 1 ) " 7ri ' 1 = ° (13) 



Except for ir , the other 7r p 's are rational functions of g and cr(g), which 
satisfies the equation [14] 18g 2 + a(l + cr)(l + 2a) = 0; in other words, 

*(9) = ~\ + ^ cos l[« ~ cos-^o) 2 ] 9l = ^ (14) 
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Starting with [14] 

" " 6s(l + 2a) 



<7(l + 3<r) 

= — 7771 777T ( 15 ) 



one finds 



ct(3 + 15ct + 16ct 2 ) a{l + a) 

772 = 2(l + 2a)3 ^ = "2(1 + 2^)3 (16) 

ct(1+ct)(1 + 5ct + 5ct 2 ) ct 2 (1+ct) 2 

773 = 3(1 + 2<t)^ ^ = (1 + 2^ (17) 

4^(1 + ^)3 

1,1,1 3(l + 2 ( r) 4 (l + 6 ( T + 6s 2 ) ff 1 ; 



and so on. 



2.4 Punctures 

So far I have been discussing how to calculate the LTE of the free energy 
F q (g, a). By definition, its singularity is F q ps t 2 ~ lsir , where t is the distance 
from the critical point (for models with multicritical behavior, this definition 
assumes that we approach the critical point from a generic direction). Now, 
since 7 str is expected to lie in the range [— |, +|], this singularity is rather 
weak, and consequently difficult to study using series methods. 

Fortunately, this problem can be circumvented by puncturing the surface. 
A puncture is a labeled, degree- 1 vertex, either in the original matrix- model 
fatgraph, or in the skeleton graph. I will denote the connected n-puncture 
Green's function G n (F q = G ). The advantage of the punctures is that 
each one brings the power of the singularity down by one: G n ps tf-n- lstr _ 
As for the matter, one has three options: putting fluctuating spins on the 
punctures, putting fixed spins, or putting no spins at all. I will take the last 
option, as it is the simplest. In this case the vertex factor for each puncture 
in the skeleton graph is just unity. 

Making three or more punctures, therefore, ensures that the critical point 
will always be a divergence. The stronger the divergence, the easier it is to 
study it using series extrapolation techniques; practically, though, making 
too many punctures is calculationally cumbersome. As a compromise I have 
chosen to make four punctures. 

Starting with the free energy, in the general case it is of course impossible 
to obtain the Green's functions. In the lte representation, however, one can 
puncture the surface in a completely mechanical way. The prescription is as 
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follows: for every partition p, replace the corresponding vertex factor n p in 
F q = G by the series 

TT p TT p + TT Ptl X + ^n PtlA X 2 H h ^^,1,...,!^ H (19) 

k 

Expand everything in powers of x: G n will then simply be the coefficient of 
x n divided by n\. The reason this works is that the topology of punctures 
is trivial; puncturing a fatgraph does not change its topology. 

2.5 Coloring the skeletons 

Each skeleton graph in the LTE contributes one last factor: the number of 
ways to assign a spin or matter state to each spin cluster, i.e., the number of 
ways to "color" its vertices. This is the only factor in the LTE that depends 
on the matter model that is coupled to the surface: the matter model defines 
the coloring rules. 

The matter model is defined by the matrix Ty that appears in the def- 
inition of the multi-matrix model, eq. (2); the entries of this matrix are 
functions of z = e~ 2t3 , the matter temperature. Ty can be considered as the 
adjacency matrix of the target space graph ("target graph" for short), T, 
each of whose vertices represents a matter state. If = then vertices i 
and j in T are not connected; otherwise they are, with the weight T i} ■,. This 
is how the coefficients of the lte depend on the target graph T: each spin 
cluster (in each skeleton) must be "colored," i.e., assigned a matter state 
(matter state = vertex of T), but only in such a way as is allowed by T and 
weighted by a correponding factor. Each skeleton contributes a sum over all 
such colorings. 

More formally: consider a skeleton graph S. A coloring of S is a mapping 
r from the vertices of S to the vertices of T. The coloring factor x(5) is a 
sum over all distinct mappings r. Each mapping is weighted by a product 
over the edges of S: an edge between vertices i and j contributes the factor 
T T (i),T(j)- Symbolically, 

X(S)= E II T rd),rU) (20) 

In other words, a mapping r is allowed only if every pair of nearest neighbors 
{i, j} in S is mapped to a pair of nearest neighbors {r(i),r(j)} in T ■ Such 
mappings are called "graph homomorphisms" or "H-colorings," and were 
introduced into graph theory by Nesetfil [15]. 
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The simplest matter model in some sense is the g-state Potts model. Its 
target graph T is the complete graph on q vertices K q , i.e., q vertices, each 
of which is connected to all the others. Every edge in T has equal weight, 
z. What are the coloring rules for this model? 

Every vertex of T is connected to every other vertex besides itself: this 
means that all mappings are allowed except ones that assign the same color to 
nearest-neighbor vertices in the skeleton. This is precisely what is meant in 
graph theory by "properly coloring" the graph S with q colors. The number 
of such proper colorings is called the chromatic polynomial of 5, x a (5): it is a 
polynomial in q of degree less than or equal to the number of vertices n in S. 
If S is a tree then the chromatic polynomial is trivially q(q— l)™ -1 ; otherwise 
it is a polynomial of degree less than n. To recapitulate: the coloring factors 
for the Potts model are just the chromatic polynomials. These are given in 
Table 1 in terms of A = q — l. 6 

I use a simple and elegant algorithm, originally due to Birkhoff, for 
computing the chromatic polynomials. Consider a skeleton graph S and 
let S' be S with some edge {«, j} deleted (provided that S' is connected). 
Deleting the edge increases the chromatic polynomial: < X g (<->') for 

any q, since there are fewer restrictions on coloring in S'. The difference is 
precisely due to those colorings for which r(i) = r(j), i.e., the colors of i 
and j coincide; in other words, if we contract i and j to a single point and 
call the resulting graph 5", we have 

Xq (S)= Xq (S')- Xq (S") (21) 

By applying the reduction (21) recursively, one eventually ends up with 
graphs that are trees, whose chromatic polynomials are trivially q(q — l)™ -1 
(n is the number of vertices of the tree). An efficient algorithm that uses 
this method has been implemented by Wilf and Nijenhuis [16]. 

Given a particular skeleton graph S and a particular target space graph 
T, one can always calculate the homomorphism coloring factor from S onto 
T. For the simple case of the Potts model (where T = K q ), we can do more: 
for a given S we can calculate, using Birkhoff 's algorithm, the coloring factor 
onto K q for arbitrary q — the chromatic polynomial of S. This means that 
the coefficients of the LTE will be explicit functions of q — one does not have 
to redo the calculation for every new value of q, and one can plug in q's that 
are not positive integers. 

6 Since my convention is to normalize matrix model expectation values by the total 
"volume" of the target graph, i.e., the number of matter states, the coloring factors in 
Table 1 are the chromatic polynomials divided by q. 
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Other models, such as multiple Ising models, are potentially more inter- 
esting (since they yield conformal theories even on the fixed lattice). The 
target graph T for v Ising models is a ^-dimensional hypercube. Its edges 
are weighted by z, corresponding to changing the state by flipping one spin; 
since more than one spin can be nipped at a time, however, the vertices 
of the hypercube are connected along the diagonals, which connections are 
weighted by higher powers of z. The varying weights in this target graph 
make the problem of calculating the coloring factor more complicated that 
the Potts case (T = K q ), as the simple Birkhoff reduction formula (21) no 
longer holds. 

2.6 Putting it all together 

This section summarizes and brings together the calculations required to 
develop the LTE. 

First, we need a list of skeleton graphs. At this stage, we consider skele- 
tons as ordinary, thin graphs. To save computational effort, the skeletons 
should be unlabeled. Initially, we make only one restriction: vertices are 
not allowed to be connected to themselves (this would contradict the def- 
inition of spin clusters); multiple edges between vertices are, on the other 
hand, allowed. 7 We will group the skeletons by their number of edges: n 
edges contributes to n-th order of the LTE. It is not trivial to generate such 
unlabeled graphs. The graphs used here have been generated by Brendan 
McKay using his software package "nauty" [17]. All the skeletons through 
order 4 are shown in the left-hand column of Table 1. Although there are 
many fewer unlabeled graphs than labeled ones, the number of unlabeled 
graphs grows faster than exponentially (there are 1183 skeletons in order 
8, 4442 in order 9, and 17,576 in order 10). Since we are using unlabeled 
skeleton graphs, each graph has a symmetry factor, the usual inverse order 
of the graph's automorphism group. These have also been calculated by 
"nauty", and are given in the "Sym." column of Table 1. 

Next, we must consider the skeletons themselves as fatgraphs. First off, 
each vertex with k legs contributes l/k\; to avoid confusion, the products of 
these vertex factorials are listed separately in the "V.F." column of Table 
1. Then recall that each vertex comes in several flavors, the twists, and 
that each twist contributes a different vertex factor: for a vertex with 3 
legs, the possibilities are 7r 3 , 7r 2 ,i, and Topologically, the twists may 

7 Graph theorists would call these "free, undirected multigraphs without slings." 
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be represented as nonlinear vertices: for 3 legs we have tr</? 3 , tup 2 tup, and 
(tr (p) z . For each vertex in a given skeleton one must sum over all the possible 
twists. For a given combination of twists (on a given skeleton graph), we 
then calculate the usual large-iV fatgraph factor: assuming that the legs 
are rigid and labeled, this is the number of ways of connecting the legs so 
that the fatgraph is spherical. For each skeleton this yields a polynomial in 
the 7r vertex factors (which are themselves non-polynomial in </), given in 
the "Spin cluster/topology" column of Table 1. The "Total" column gives 
the product of the symmetry, vertex factorial, and spin cluster /topology 
factors, expressed using the convenient p vertex factors (which are just linear 
combinations of the 7r's — see section 2.2). Finally, in order to ease the 
subsequent numerical analysis, the surfaces need to be punctured. This is 
done using eq. (19) — see section 2.4. In the analysis in the following section, 
I will make four punctures. 

In order to specialize to a particular target space (matter model), one 
must calculate a homomorphism coloring factor for each skeleton onto the 
target space — see section 2.5. For the g-state Potts model, this coloring 
factor is just the ordinary chromatic polynomial. Through order 4, these 
are given in the right-hand column of Table 1, normalized by q and expressed 
in terms of A = q — 1. (Recall that a = Az.) 

In the present work I have carried out the LTE for the Potts model 
through order a 10 . Computationally this has not been very taxing. Using 
rather elementary algorithms implemented in Fortran, the whole calculation 
took about two hours of computer time on an HP 9000/700 workstation. 
(Actually plugging in various values of g and q into the expressions obtained 
for the coefficients also takes a non-negligible amount of time.) 

The resulting series have been checked in two ways. One can set q = 2 
or q = oo and compare to a low-temperature expansion of the two available 
exact solutions. Alternatively, one can carry out, by hand, an expansion in 
powers of g for any q (the so-called strong-coupling expansion) up to order 
g m (the coefficients will be polynomials in a); grouping together powers of 
a, we obtain an approximate LTE to order « 3m / 2 (approximate because each 
coefficient — a non-polynomial function of g — is given as a Taylor expansion 
to order g m ). 
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3 Example: the Potts model 



3.1 Preliminaries 

Series extrapolation is always a tricky business, and especially in the present 
case, where each coefficient of a n is itself a complicated function of g and 
q. Fortunately, there are several techniques available (the Pade and ratio 
methods, and countless variations thereon), which allows for cross-checking; 
all results to be presented are supported by more than one technique. In 
addition, much is known exactly for the cases q < 4 and q = oo, which also 
provides a valuable check. 

Near its critical point a c (q,g), I will assume that the 4-puncture function 
has the form 

G 4 (q,g,a) « \a - a c (q,g)\- 2 -^^(\og \a - a c (q,g)\) x ^ (22) 

(up to multiplicative and additive analytical corrections). The only reason 
to assume this form is that it is correct for the exactly solved cases (q < 4 
and q = oo). The form of the logarithmic correction is really a guess; it is 
exact for the d = 1 model, where j str = and A = — 1. Expanding in powers 
of a 

G 4 (q,g,a)=Y^c n (q,g)a n , (23) 

n 

one finds asymptotically for high orders 

c n (q, g) ~ a c (q, g )- n 1+ ^^ (log n) x ^ (24) 

For a given q, there is a critical curve (g,a c (q,g)), on which the sur- 
faces become large (the continuum limit); on this curve there might be one 
(or more) multicritical points (g*(q), a*(<?)) which correspond to the phase 
transition(s) of the matter. A quantity that will be useful is the free energy 
density of the matter in the limit of large surfaces. This can be obtained as 
follows: first rewrite (23) as a series in g: 

Gi(q, g,a) =J2 d n(<l, a)g n (25) 

n 

with the asymptotics of the coefficients 

d n (q,a) ~ gM^r n n 1+lsU9 ' a \^n) K9 ' a) (26) 

The n-th term of the strong-coupling expansion (25) counts surfaces of area 
n: the partition function of the matter in the limit of large area (n — > oo) 
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is then precisely l/g c (q, a), which can be obtained by inverting the function 
a c (q,g), calculated using (24). The free energy of the matter, 



U(q,(3) 



log g c (q, a), (3 



log a/ (q - 1) 



(27) 



should have the following scaling form near a multicritical point (g*,f3*) 



It is important to distinguish the G 4 scaling relation (22) from the U scaling 
relation (28). In the first we approach the critical (</, a) curve from some 
generic direction; while in the second we are constrained to that curve, 
and approach the multicritical point (</„,, a* ) . The U scaling may also have 
logarithmic corrections (as for the q — > 1 case, where U ~ (/3„ — (3) 2 log(/3* — 



The crucial problem is locating the phase transition. First of all, for 
large q, there are probably two transitions rather than one: between the 
low-temperature magnetized phase and the branched polymer phase; and 
between the branched polymer phase and the high-temperature disordered 
phase. The high-temperature transition seems rather hard to detect, how- 
ever. Both in the lte and in Monte Carlo experiments [18], the transition 
from the branched polymer phase to the disordered phase seems gradual, 
though there are some very good theoretical reasons that it not be grad- 
ual [19]. Most likely, the finite-size effects smooth out the high-temperature 
transition much more than the low-temperature one; though why this should 
be so is a mystery. In any case, the phase transition that I will discuss in 
the rest of this paper will be the low-temperature one: the one out of the 
magnetized phase, whether it be into the disordered phase (for low q), or 
into the branched polymer phase (high q). 

I will use two methods to locate the phase transition. First, for q < 
4, I will assume that there are only two phases, the magnetized and the 
disordered (both with ^y str = — |), and at the phase transition ^y str > —\- So 
the way to locate the phase transition, as well as to calculate 7 sir , is to find 
the maximum of j str along the critical curve (the continuum limit). This 
method is used in conjunction with Pade approximants — see section 3.2. 

For larger q this method is unreliable, however, due to the probable 
presence of the branched polymer phase (which should have 7 sir = +|) 
between the ordered and the disordered phases (which still have j str = 
— |). Therefore the exponent at the phase transition — | < ^y str < +| is 
no longer the maximum j s t r - Another way to find the phase transition 



t/(g,/3)~|/3-&| 2 -^ 



(28) 



£))■ 
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is to use the peaks of the of the derivatives of J7, U k = d k U/d(3 k ; this 
is used in conjunction with the ratio method — see section 3.3. The peak 
in the specific heat, U 2 — and in all higher even derivatives, U 4 ,U 6 ,... — 
should indicate the phase transition, with one exception. Assume that the 
exponent a is an integer, n, and that there are no logarithmic corrections to 
the scaling relation, eq. (28). This is indeed the case for q = 2 and q = oo, 
where a = —1. This means that U, Ui, . . . , ?7i_„ are continuous at the 
phase transition, and the slope on either side is finite; the higher derivatives 
U-i-n-, ■ ■ ■ are discontinuous at the critical point, but again with finite slopes. 
Take a = —1, for example: although the specific heat U 2 might still have a 
peak, this is not necessarily true, and will not be the case if the slopes on 
either side of the critical point have the same sign. The next derivative, U 3 , 
will have a discontinuity, but since there are no divergences this does not 
have to result in a peak, or it might give two peaks, one on either side of 
the transition. It is only U 4 that must have a peak roughly at the critical 
point. In general, for a = n with no logarithmic corrections, ?7 3 _„, £7s_ n , . . . 
should have single peaks near the critical point, while U 2 -m U 4 - n , . . . should 
have double peaks, one on either side of the critical point. 

To summarize: if a is not an integer, or if it is an integer but with 
logarithmic corrections, then U 2 ,U 4 ,... should have peaks roughly at the 
critical point. If a is some integer n but without logarithmic corrections, 
we should look for peaks in U 3 _ n , ?7 5 _„, ... to indicate the phase transition. 
What is the likely range of a? Between q = 2 and q = 00 I will assume that 
a > — 1, the endpoint values. We know from KPZ scaling that a(q = 4) = 
(it is not known whether there should be logarithmic corrections). Prom 
a Monte Carlo study [18] we know that at a(q = 10) is very close to 0, 
while at q = 200 it is down to roughly —0.85. It seems that in all plausible 
cases, therefore, U 4 should have a peak near the critical point, except ifa = 
without logarithmic corrections, in which case the U 3 peak should be a better 
indicator of the critical point. Not wishing to make a priori assumptions 
about a(q), I will therefore present two candidates for the critical point, the 
locations of the U 3 peak and the U4 peak. It will be seen that the U4 peak 
fits well for small and large q, while the U 3 peak makes more sense for an 
intermediate range of q > 4, as expected. 

3.2 Pade approximations 

The Pade method [20] assumes that A = 0. If there are, indeed, logarithmic 
corrections, what one obtains for 7 str is an effective exponent which is shifted 
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from the true value in the same direction as the sign of A. One fits the 
logarithmic derivative of (22) 



|- log G^g.fra)* , 2 7 7 (29) 
oa \a — a c \ 

by a rational function, degree-L in the numerator and degree-M in the 
denominator 

[L, M] = Po+P^ + Pif + ---+PL* L m 

1 +qia + q 2 a 2 + ••• + q M a M V ' 

The coefficients can be calculated by simply expanding the right-hand side of 
(30) in powers of a, and equating to the logarithmic derivative of G 4 (though 
there are more effecient methods). Therefore, if the series is known through 
order a n , once can calculate all the [L, M] provided that L + M < n — 1. The 
zeros of the denominator of (30) give the locations of the singularities (some 
of which are physical and some of which are spurious), and their residues 
give the critical exponents. By constructing a number of different [L, M] 
approximants, one can get a sense for which are the physical poles (the ones 
that recur), and a subjective impression of the error. As is often the case, 
the [L, L — 1] and [L, L] approximants seem to be most accurate. 

As an example of the Pade method, Fig. 4 shows the [5, 4] approximant 
for the exponent 7 sir for the physical pole as a function of the cosmological 
constant g, for q = 2, 4, 10, 25, 100, 500, and 10000 (as q gets larger, the 
peaks get higher). The right-hand side of the graph corresponds to low 
temperature and pure gravity: the magnetized phase, hence ^y str = — \. The 
curves were cut off on the left when a c reaches A = q — 1, i.e., when the 
matter temperature reaches infinity. For a given value of q, as one varies g 
the locations of the physical and spurious poles vary. When a spurious pole 
crosses the physical pole, one sees a little "glitch" : a localized divergence in 
7 sir as a function of g. In Fig. 4, such glitches have been smoothed over to 
make the plots easier to read. 

For low q, certainly for q < 4, we can interpret the curves as follows: the 
flat, 7 sir pa — I regions on the right- and left-hand sides of the peak corre- 
spond to the pure gravity, magnetized and disordered regimes, respectively: 
while the peak corresponds to the phase transition, at which the divergent 
correlation length of the spins (for q < 4, at least) modifies the global ge- 
ometry of the surface and therefore j str (or alternatively, the point at which 
the flat conformal field theory is dressed by gravity). According to the KPZ 
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Figure 4: [5,4] Pade estimate of ^y str as a function of the cosmologi- 
cal constant </, for different values of q (in order of ascending peaks, 
q = 2,4, 10,25, 100,500, 10000) 
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Exact 


[4,3] Pade 


[5,4] Pade 


9c 


a c 


Ifstr 


9c 


a c 




9c 


a c 




2 


0.0694 


0.159 


-1/3 


0.0703 


0.135 


-0.371 


0.0699 


0.146 


-0.358 


3 


0.0677 


0.254 


-1/5 


0.0687 


0.222 


-0.286 


0.0682 


0.240 


-0.266 


4 









0.0672 


0.303 


-0.204 


0.0673 


0.300 


-0.207 



Table 2: Results from Pade analysis for low q using j str peaks to locate the 
phase transition 

formula, valid for q < 4, the value of j str at the transition is 

The results for low values of q are presented in Table 2 (the exact a c and 
g c are taken from [21, 22]). The results are quite different for q = 2 and 
3 than for q = 4. While the location of the critical point and 7 sir are 
not very accurate for q = 2 and 3, they are clearly moving in the right 
direction as one goes from the 8th order of the lte (the [4,3] Pade) to the 
10th order (the [5,4] ) 8 . For q = 4, on the other hand, the results hardly 
change at all from 8th to 10th order; j str actually gets a little worse. This 
"sluggishness" is symptomatic of the following general phenomenon: when 
there logarithmic corrections to scaling are present the convergence of Pade 
(and other) approximants is drastically slowed. It should also be noted that 
a value of 7 sir ps —0.2 is commonly obtained in Monte Carlo simulations for 
various c = 1 models. 

In Fig. 4 one clearly sees a broadening of the 7 sir peak, as has been 
observed in Monte Carlo experiments [23] and in strong-coupling series study 
[7]. For the larger q, what we have no longer looks like a peak, but rather 
a plateau. Could this be the branched polymer phase, predicted in [9] to 
lie between the low-temperature (magnetized) and the high-temperature 
(disordered) phases? The presence of a plateau (which extends all the way 
to g = for q = oo) for large q means that the "peak" in j str can no longer 
be used as a way to detect the phase transition. Indeed, the "peak" values 
of 7 s(r clearly approach +| as q — > oo (the branched polymer value), while 

8 Because of the logarithmic derivative (29), given the series to order a n (including the 
zeroth order term), we can calculate Pade approximants up to L + M = n — 1 
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we know that at the phase transition, j str = | when q is sufficiently large. 
In what follows, I will use peaks in the specific heat and its derivatives to 
locate the phase transition. 

Although it gives the correct general picture, there are three reasons why 
the Pade method is ultimately not very useful. First, it is difficult to modify 
it to include logarithmic corrections to scaling (A / 0). Second, there is 
no general procedure to extrapolate from finite [L, M] to L, M — > oo. And 
finally, as the glitches due to root crossings have a strong and random effect 
on a c (and a stronger effect on its derivatives), one cannot reliably use the 
Pade method to look for peaks in the specific heat (and its derivatives) in 
order to locate the phase transition for high q. 

3.3 Ratio method 

A more adequate procedure in this case is the ratio method [20]. The first 
step is to take the ratios of the coefficients c n (q,g) in (23), which should 
have the asymptotics (24): 

r n = -5- ~ - ( n \ 1+1StT ' ( l ° gn ) X (32) 
c„_i a c \n-l) Vlog(n-l)/ 

Typically, one expands the right-hand side of (32) in inverse powers of n, 
the constant term giving a c , the coefficient of 1/n giving 7 s(r , and the coeffi- 
cient of l/(n log n) giving A; these can then be extrapolated to n — > oo, their 
extrapolants further extrapolated, etc., giving successively more refined esti- 
mates. This procedure works for the LTE, but I will adopt a slightly different 
method which seems to be more stable and, for the cases where the answers 
are known, more accurate. 

For each value of n (and q and g), calculate the "running" values of 
a c (n;q,g), 7 str (n; q, g), and X(n;q, g) by solving eq. (32) for n, n — 1, and 
n — 2. Assuming the scaling form (22), with multiplicative and additive 
corrections that are analytic at a c (q,g), the running values should approach 
their n — > oo limits with corrections that are 0(1/ n); and, if logarithms 
are present, with additional corrections of C(l/nlogn). Fitting the running 
values to the above form, one estimates the asymptotic values. 

The are two main questions to be answered: are there logarithmic cor- 
rections to scaling, and where exactly is the critical point? If logarithms are 
indeed present, we may safely fit to the asymptotic form (32). If not, how- 
ever, fitting to logarithms can be quite dangerous. We know that there are 
no logarithms for q < 4 and q = oo; and indeed, if one fits to (32) there, one 
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Figure 5: Values of ^y str obtained as a function of q, using different methods 
to locate the critical point, with and without logarithmic scaling (see text). 
The KPZ values are also shown for q < 4. 



gets nonsense. So the question is, for which values in the range 4 < q < oo 
should one fit to logarithms? In what follows, I will hedge my bets, fitting 
to (32), as well as to the same form with A = 0. I will locate the critical 
point using peaks in U 3 and U 4 . As discussed above, the U 3 peak should 
be used only if a = without logarithmic corrections, and otherwise the U 4 
peak should be used. I will present data at both peaks. 

Fig. 5 presents values of ^y str calculated for a wide range of q using the 
methods discussed above. The phase transition is located using the peak in 
U 3 or U 4 ; and the data are fitted with and without logarithmic corrections, 
For some values of q, more complete data are given in Table 3. 

There seem to be three different regimes: low q (q < 4), intermediate q 
(4 < q < 50, rougly), and high q (50 < q < oo). First, low q. According 
to KPZ, a = — 1 (discontinuity) for q = 2 and a = — \ for q = 3 (so we 
should use the peaks in U 4 for locating these transitions), and there are no 
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0.131 


-0.347 


0.159 
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0.206 


-0.267 


0.251 


-0.299 


4 


0.290 


0.00329 


-0.710 


0.352 


0.0384 


-0.780 


0.260 


-0.206 


0.319 


-0.248 


5 


0.337 


0.0359 


-0.711 


0.408 


0.0219 


-0.724 


0.301 


-0.158 


0.370 


-0.204 


6 


0.374 


0.0741 


-0.726 


0.454 


0.0488 


-0.725 


0.334 


-0.121 


0.413 


-0.168 


10 


0.476 


0.169 


-0.753 


0.585 


0.104 


-0.706 


0.423 


-0.0231 


0.530 


-0.0679 


25 


0.627 


0.200 


-0.644 


0.796 


-0.0278 


-0.421 


0.553 


0.107 


0.714 


0.0934 


50 


0.700 


0.102 


-0.427 


0.916 


-0.304 


-0.0345 


0.614 


0.158 


0.812 


0.196 


100 


0.737 


-0.0294 


-0.194 


1.018 


-0.736 


0.493 


0.644 


0.174 


0.870 


0.268 
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0.665 


0.166 


0.925 


0.360 


oo 














0.666 


0.162 


0.930 


0.372 



Table 3: Critical points and exponents calculated by the ratio method 



logarithms (i.e., A = 0). For q = 2, the U 4 peak gives a very good value of 
a c = 0.1586, compared to the exact answer, a c = 0.1589, as well as a decent 
value for ^y str . For q = 3, a c is still very close to the true value a c = 0.254; 
7 sir , on the other hand, is not very good. Part of the problem is that the 
peaks in j str are quite narrow (which is a sign that the method is working 
well), so a small error in the location of the critical point results in a large 
error in 7 str . This problem should be less acute for larger q, where the peaks 
are broader. In any case, the method gives quite convincing results for the 
location of the critical point, which is what we want from it. 

Although it is not known exactly, one suspects that for q = 4 (c = 1) 
and higher, there are logarithmic corrections to scaling: that is, A ^ in 
eq. (22). Therefore I do the fit (32) with the logarithm. We have a = at 
q = 4 by KPZ, and a ?a at q = 10 by Monte Carlo; but we do not know 
whether a is also modified by logarithms or not. If it is, then following the 
discussion in section 3.1 one should use the peak in J7 4 to locate the phase 
transition; but if a = is not modified by logarithms, one should use the 
U 3 peak instead. Consider the q > 4 region of Fig. 5. A striking result 
is the difference between the logarithmic and the non- logarithmic fits. At 
q = 4, both non- logarithmic fits give j str « —0.2, in agreement with Monte 
Carlo measurements, and in disagreement with KPZ. Just by including a 
logarithmic term one obtains radically better agreement with 7 sir = 0. Since 
the U 3 j s tr is an order of magnitude closer to zero than the U 4 j s t r , we have 
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evidence that at q = 4, a = without logarithmic corrections. Thus for 
the c = l,q = 4 model I conclude that ^y str pa in agreement with KPZ, 
but with logarithmic corrections: A pa —0.7. This should be compared to a 
Monte Carlo calculation for the same model which gives A pa —1.5, as well 
as another c = 1 model, the continuous d = 1 model, which has A = — 1. 

For q > 4 the logarithmic corrections are likely to persist. Up until 
q « 10, all four estimates for ^y str grow. The U 3 and U 4 non- logarithmic 
estimates follow each other closely, and cross zero around q = 10 — 15. The 
logarithmic approximations also increase, the U 3 faster than U 4 . However, 
by q pa 15 (for U 4 ) and q pa 25 (for U 3 ) the 7 sir turns around, eventually 
reaching zero and becoming negative. The logarithmic exponent A stays 
negative and rather flat in the range q pa 4 — 10, and then goes to zero, 
around where ^y str begins to decrease. For the U 3 case, A seems to go to 
zero, while for the U 4 case it becomes positive. As I have already mentioned, 
logarithmic fits are dangerous when there is no logarithmic scaling, or when 
the logarithms are weak. This is precisely what happens here: when the 
logarithms start to weaken, the entire approximation starts to get erratic. 
It seems that in the range q = 4 — 15 or 20 there is a rather flat logarithmic 
correction, with the exponent A pa —0.7 or —0.8; after which the logarithmic 
correction decreases (it is impossible to say at present whether this decrease 
is gradual or sudden), and one can no longer trust the logarithmic fits. 

Around q pa 20 for the U 4 case, and q pa 40 for the U 3 case, the log- 
arithmic estimates for j str cross their non- logarithmic counterparts. Since 
all indications are that the logarithmic corrections fade out by that point 
(remember, at q = oo we have A = 0), we should trust the non- logarithmic 
estimates above q pa 40. At this point the non- logarithmic U 3 estimate di- 
verges from the U 4 estimate. As q — > oo, we have ^y s t r (U 3 ) — > 0.16, while 
lstr{U 4 ) — > 0.37. Only the latter fits well with the exact solution of the 
q = oo model [8], where a = — 1 (without logarithms: therefore the peak in 
U 4 should signal phase transition) and 7 sir = 1/3 at the phase transition. 
All of this indicates that above q pa 40 the non- logarithmic U 4 estimate for 
the critical point is the one to use. 

It might seem worrisome that 7 sir overshoots its q = oo value 1/3, 
settling down to the somewhat higher value 0.37. Similaly, a c — > 0.93 
(in place of the exact a c = 1), and g c — > 0.0582 (in place of the exact 
l/\/288 pa 0.05893). This is, however, a finite-order artefact of the LTE. 
This can be checked by expanding the exact q = oo solution in LTE. What 
happens is that the approximants are slightly oscillatory (with period two: 
a common phenomenon, usually due to nearby unphysical singularities, and 
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which is often cured by an Euler transform in the complex a-plane; but the 
cure can sometimes be worse than the disease), and don't completely settle 
down until n > 10. An ad hoc solution is to take the alternating approxi- 
mants a(8) and a(10) and extrapolate them linearly in 1/n to n — > oo, and 
thereby obtain a more refined estimate for U and its derivatives. This gives 
a much better location for the critical point: a c = 1.011, g c = 0.05898. 

If the magnetization phase transition were first-order (as it is on a flat 
surface for q > 4), the exponent j str at the muticritical point should be 
exactly — |, the pure gravity value: only critical fluctuations should modify 
j s tr- Except possibly for the strange region around q = 40, j str gives no hint 
of descending back down to — |, for any value of q. One therefore concludes 
that coupling a Potts model to a random surface softens the phase transition, 
making it continuous for most — if not all — values of q. 9 

4 Prospects 

In this paper I have shown how to develop the LTE for any multi-matrix 
model on any target space. Computationally, the task is to calculate ho- 
momorphism coloring factors onto the desired target space. For the case of 
the Potts model, these homomorphism coloring factors are just the ordinary 
chromatic polynomials which are very easy to calculate. So for this model 
I have worked out the lte out to 10th order. This takes about two hours 
of computer time on a workstation, and gives quite accurate results for not- 
so-branched (q < 4, c < 1) surfaces, as well as for very branched (large q, c) 
ones; and a plausible account of the mysterious region in between. 

It is important to find ways to carry the LTE to higher order than 10th 
(the fact that it takes about 2 hours of computer time to go this far is 
promising in this respect). This is not just busywork. In ordinary spin 
models, especially with logarithmic corrections, one starts to get very accu- 
rate estimates for the exponents in order 15 or 20, typically. A similar effect 
can be seen in the q = oo model, which, since it can be solved exactly, can 
be expanded to arbitrary order in the LTE: the series is rather chaotic until 
order 6 or 7, then starts to behave well, but does not become completely 
smooth until order 15 or so. Thus it could be very profitable to push the 
LTE a few orders further. One way to do it would be to count only 1- or 
2-particle irreducible skeletons; the full series can then be reconstructed. 

9 A possibly interesting comparison can be made to the Potts spin glass, which, for 
sufficiently large bond disorder, also has continuous transitions for all q [28]. 
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For example, of the 17,576 skeletons in 10th order, only about 200 are 2PI. 
Another improvement would be to calculate moments of the magnetization 
(rather than just closed surfaces or punctured ones with no spins on the 
punctures). One could then use the powerful Binder's cumulant technique 
to calculate the location of the critical point, rather than the peaks of the 
specific heat and its derivatives. 

Having shown that the lte works well for the Potts model on a random 
surface, nothing stands in the way of extending the calculations to other, 
more interesting, models. Well, almost nothing: one first has to devise a 
way of computing homomorphism coloring factors onto general target space 
graphs. But this is only a computational, rather than a conceptual, problem. 
Examples of models to study are d-dimensional lattices and multiple Ising 
models, which can easily implement c > 1 (for some recent work, see [24, 
25, 23, 26, 27]. At c = oo these models are identical to the q = oo model 
[4], but the region of intermediate central charge is still poorly understood. 
The LTE seems to be a promising tool for its investigation. 

I am grateful to LPTHE, Jussieu for its accueil chalereux and to Marco 
Picco for useful conversations of all kinds. The graphs used in this work 
were kindly provided by Brendan McKay. 
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